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Anti-solar differential rotation 
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Abstract. The differential rotation of anti-solar type detected by observations for several stars may result from a fast merid- 
ional flow. The sufficiently intensive meridional circulation may be caused by large-scale thermal inhomogeneities or, per- 
haps, by tidal forcing from a companion star. First results of simulations of the anti-solar rotation of a giant star with mag- 
netically induced thermal inhomogeneities are presented. Perspectives for observational check of the theory are discussed. 
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1. Introduction 

This paper discusses possible theoretical explanations for ob- 
servations of non-uniform stellar rotation with the angular 
velocity increasing with latitude. It also suggests the possi- 
bilities for observational check of the theory. 

The anti-solar rotation is relatively seldom to observe. 
The solar-like case of rotation rate decreasing with latitude 
is more frequent to detect (Strassmeier 2003 Petit et al. 



2004). The relatively fast rotation of the equator represents 
also the case which theory succeeded to explain. Theory of 
angular momentum transport by rotating turbulence and nu- 
merical simulations both show meridional fluxes of the an- 
gular momenrtum towards equator (Riidiger & Hollerbach 



2004). Models for differential rotation based on that the- 
ory provide rotation law for the Sun (Kitchatinov & Riidiger 
1995, Kitchatinov 2004) in close agreement with helioseis- 
mology (Shou et al. I1998> . Observations of solar-type ro- 
tation for AB Dor (Donati & Collier Cameron [19971 and 
LQ Hya (Kovari et al. 2004) can be reproduced with reason- 
able accuracy (Collier Cameron et al. 2001 Kitchatinov et al. 
2000), and theoretical trend for differential rotation of main- 
sequence dwarfs to decrease with spectral type (Kitchatinov 
& Riidiger 1999 ) was given some observational support (Col- 
lier Cameron l200T1 Petit et al. 120041 ). 

However, observations of anti-solar rotation are already 
numerous enough to demand a revision of theory. The stan- 
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dard hydrodynamical models cannot reproduce this case. The 
most plausible way for such a revision is to include some ad- 
ditional driver of meridional flow. We shall see in the next 
section that the fast meridional circulation is the most clear 
theoretical possibility for producing differential rotation of 
the anti-solar type. Observations may help a lot in under- 
standing the nature of stellar differential rotation by checking 
whether all the stars which show anti-solar rotation do indeed 
possess a fast meridional flow. 

The fast meridional motion can result from a barocline 
driving due to large thermal spots or from tidal forcing by 
a close companion. The guess is suggested by both theoret- 
ical arguments and statistics of differential rotation observa- 
tions for individual stars summarized by Strassmeier (2003). 
All stars with detected anti-solar rotation belong to one of 
two (yet small) groups: (i) close binaries, or (ii) rapidly ro- 
tating giants with large-scale thermal inhomogeneties (Hack- 
man et al. 2001 Strassmeier et al. 2003 ). The suggestion that 
all anti-solar rotators are spotted giants or close binaries is 
more speculative and less certain than the conclusion on the 
role of meridional flow but may also be a subject for obser- 
vational verification. 

After a general discussion of possible origin of anti-solar 
rotation in the next section, section [3] outlines a numerical 
model for equatorial deceleration on a giant star with large- 
scale thermal inhomogeneties caused by magnetic fields. Re- 
sults of first simulations are presented and discussed in sec- 
tion|4] The final section[5]suggests the ways for observational 
verification of the theory. 
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2.3 Convective heat transport in magnetized fluids 



2. General consideration 

2.1. Anti-solar rotation for fast meridional flow 

Possible reason for anti-solar differential rotation can be in- 
ferred from consideration of a global axisymmetric flow in a 
spherical convective shell. The flow velocity, u, can be ex- 
pressed in terms of the stream-function, ip, of the meridional 
circulation and angular velocity, ft, using the standard spher- 
ical coordinates, r, ip, 

1 dip — 1 dip 



u = 



• sin i? Q 



(1) 



v pr 2 sin i? d'd ' pr sin $ dr ' 

The flow obeys the steady motion equation, 

p (u ■ V) u + ((B • V) B - VB 2 /2) I (4tt) 

+ Vp - pV$ = -div OQ) . (2) 

where B is magnetic field, $ is gravity potential, and Q is 
the correlation tensor of fluctuating velocities, u', 

Qij = ( u i u j)- (3) 
The most promising for producing anti-solar rotation is a 
fast meridional flow. One can find that from the equation for 
angular velocity which results as (/(-component of Eq. 

1 dip d(sm 2 m) 1 dtp d(r 2 n) 



sin 2 d dr d'd r 2 di) dr 

where only the contribution of meridional flow is written ex- 
plicitly. The dotted terms in @ can be neglected if the flow 
is fast enough. Then, the solution can be found, 



0. 



(4) 
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which describes almost certainly the inhomogeneous rota- 
tion of the anti-solar type. The angular momentum is con- 
served along stream lines of a fast meridional flow (Rudiger 
1989). F in equation (01 is an arbitrary function. Additional 
conditions are required to define it. Whatever the conditions 
could be, the solution (0 tends to describe the anti-solar ro- 
tation. Indeed, the stream-function, ip, is constant along the 
stream lines. Therefore, the angular velocity (0 varies along 
the stream lines in such a way that it increases when the lines 
approach the rotation axis at high latitudes. The conclusion 
does not depend on whether the flow is poleward or equa- 
torward (on the surface), in any case the fast flow tends to 
produce the anti-solar differential rotation. 

2.2. Barocline flow 

A barocline meridional flow can be sufficiently fast for sup- 
porting the anti-solar rotation. The barocline term contributes 
the equation for a steady meridional flow (Kitchatinov & 
Rudiger fT999l . 



2>(V0 = ^(VpxVpl 



(6) 



when the surfaces of constant density and pressure do not co- 
incide. In equation (0, the left part accounts for resistance to 
meridional circulation by eddy viscosity and only the baro- 
cline source is written explicitly in the right. It is convenient 
for our purposes to express it in terms of specific entropy, 



S, and gravity potential, as it was done by Kitchatinov & 
Rudiger ( fT99^ . 
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The entropy and gravity distributions are normally close to 
spherical symmetry and the resulting barocline flow is small. 
This is why the advection dominated states did not emerge 
and anti-solar rotation was not found in former simulations. 

Meridional circulation can be fast, however, if consid- 
erable deviations from spherical symmetry are available in 
gravity or temperature distributions. The asymmetric grav- 
ity is present in binary systems and asymmetric temperature 
is typical of giants with their large thermal spots. Anti-solar 
rotation can be expected for these cases. The guess agree 
with statistics of anti-solar rotation detections summarized 
by Strassmeier (2003 1. Between nine anti-solar rotators, six 
belong to close binaries and two - to giant stars with large 
thermal inhomogeneities on their surfaces, the remainder star 
is LQ Hya for which different observations disagree on the 
sense of its differential rotation. 

Later on, we focus on the case of thermal asymmetry to 
treat it in a more quantitative way. 

2.3. Convective heat transport in magnetized fluids 

Stellar spots are believed to be magnetic by origin. To ac- 
count for the magnetic field influence on thermodynamics, 
we include magnetic quenching of convective heat flux, 

F?™ = -pT Xij —, Xij = ^(P)x^ 



(8) 



where x° is the eddy thermal conductivity for nonmagnetic 
case and ip is the quenching function of the normalized field 

strength f3 — B / \J A-Kp(u' 2 ). The function ip is steadily de- 
creasing with (3. The quenching function was derived in the 
paper by Kitchatinov et al. ( 1994 1 where it is given as p> x - 
f unction. 

3. The model 

Our present model is very close to its previous version 
(Kitchatinov & Rudiger 1999). We describe the model only 
briefly focusing on where it is different from the former for- 
mulation. 

The main difference is that we include now axisymmetric 
magnetic field governed by the steady induction equation, 



V x (u x B - r/ T {(3) V x B) = 0. 



(9) 



Eddy magnetic diffusivity is defined in terms of superadia- 
baticity of the stratification and includes quenching by mag- 
netic field and rotation, 
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where I and r are mixing length and time respectively, the 



function <p is given in Kitchatinov et al. ( 1994i, the magnetic 
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quenching function, tp is the same as before, and Q* is the 
Coriolis number, 

n* = 2rn. (ii) 

Anisotropy of magnetic diffusion is neglected in (|9j. We 
neglect also the anisotropy of eddy viscosity, but keep the 
rotationally induced anisotropy of thermal conductivity be- 
cause the differential rotation models cannot function nor- 
mally without it (Riidiger et al. 2004 1. Magnetic quenching of 
all the turbulent transport coefficients including the A-effect 
(Riidiger 1989 1 was described by same quenching function 
f(0) ©■ The axisymmetric magnetic field can be expressed 
in terms of the toroidal field B and potential A of the poloidal 
field: 

1 OA 1 OA (12) 



Table 1. Parameters of the star 



B 



,B 



r 2 sin ■d eW r sin Or ' 

Our numerical model solved a system of five joint equa- 
tions for angular velocity, meridional flow, entropy, poloidal 
and toroidal components of the magnetic field. 

The induction equation l|9) does not include the a-effect 
(Rrause & Radler 1980) of turbulent dynamo. Therefore, our 
model cannot support any dynamo. The magnetic field was 
involved through the boundary condition of a steady radial 
field penetrating the convection zone at the inner boundary, 
Ti, from the radiative core. The bottom field was prescribed 
by the steady potential A at the inner (bottom) boundary, 

<j> m 

'S 2 "tf), Bo = — , (13) 



.4 



1 



2nrf 



where <£> m is the magnetic flux of the dipolar field per hemi- 
sphere, n is the parameter controlling the latitudinal distribu- 
tion of the field, the larger is n there more concentrated to 
poles is the poloidal field. Computations for various n were 
performed. Unless otherwise stated, the results of the next 
section correspond to n = 3. The condition Jl 3i for the in- 
ner boundary can be understood as penetration of a relic field 
stored in the radiative core into the convection zone. In con- 
vection zone proper, the field is subject to turbulent diffusion 
and advection, so that toroidal field is produced by the differ- 
ential rotation. 

The other boundary conditions were zero radial velocity 
and zero stress at the top and bottom, constant heat flux at the 
bottom and black-body radiation at the top, superconducting 
condition at the bottom for the magnetic field and vacuum 
condition on the top. 

4. Results and discussion 

We performed our simulations for a giant star evolved from 
the main-sequence. The parameters of the star were taken 
form an appropriate model for stellar structure (Herwig et al. 



1997), some of them are given in the table. The assumed ro- 
tation rate is marginal for Doppler imaging. The rate is, how- 
ever, several times larger than normal for this type of stars 
(Grav fT989l 

The differential rotation simulated for the nonmagnetic 
state is shown in Fig.^ The standard equatorial acceleration 
of about 10% was found for this case. 
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Fig. 1. Differential rotation for nonmagnetic case is of the so- 
lar type. Angular velocity increases from poles to equator by 
about 10%. 



Figure |2] shows how the global surface flow varies with 
the amplitude of prescribed poloidal field Jl 3i for several lat- 
itudinal profiles of the field. As the field increases, the merid- 
ional flow reverses to poleward orientation and then grows 
steadily. When the flow becomes sufficiently fast, the dif- 
ferential rotation reverses to anti-solar case. We were able 
to follow the dependencies up to the value of order Bq ~ 
100. ..300 G depending on the bottom profile dl3t of the 
poloidal field. Our numerical code based on the relaxation 
method did not converge for stronger fields. The numeri- 
cal instability, probably, signals on the onset of the physical 
flux-concentration instability (Kitchatinov & Mazur 120001 
via which the cool magnetic spots are formed. By this reason, 
we did not have true magnetic spots but smooth magnetically 
induced thermal inhomogeneities in our simulations. 
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Fig. 2. Left panel: surface differential rotation, 
100 (Q cq - fiP°^| /n , as function of the poloidal field 
amplitude, Bq J 1 31 . Actually computed points are marked 
by signs. Different signs correspond to different values of 
the n-parameter of the bottom profile M2>\ of the poloidal 
field: n = 1 (diamonds), n = 3 (crosses), n — 5 (stars). 
Right panel: Surface meridional flow at 45°-latitude. Positive 
values mean poleward flow. 



Figure [3] illustrates typical case of anti-solar rotation of 
present simulations. Angular velocity increases from equa- 
tor to pole but not steadily. In agreement with above qualita- 
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tive arguments, the increase occurs in the same latitude range 
where relatively fast poleward meridional flow can be ob- 
served in Fig.0] The surface rotation of Fig.[3]cannot be ap- 
proximated by traditional cos 2 i? profile. It may be reasonable 
to use higher-order terms in approximations of the observed 
anti-solar rotation laws. Expansion in terms of Gegenbauer 
polynomials (Rudiger 1989 1, 



JV 



^) = E-V2„-l(cOstf), 

^— ' smv 



(14) 



may be convenient, especially when rotation law variations 
with time are concerned. 
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Fig. 5. Left panel: lines of the poloidal magnetic field for 
same model as Fig. [3] Right panel: surface field strength. 
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Fig. 3. Simulated differential rotation for poloidal field ampli- 
tude Bo — 200 G. Angular velocity increases from equator 
to pole by about 10%. Low and high latitudes have opposite 
sense of the latitudinal rotation inhomogeneity. 



The meridional flow of Fig. |4] is driven by the barocline 
forcing due to the thermal depletion which occupies same re- 
gion of latitudes. As was already mentioned, there were no 
true thermal spots in the simulations. Nevertheless, the ther- 
mal depletion of Fig. |4] was produced by magnetic quenching 
of convective heat flux of Eq. l|8). 





Fig. 4. Left panel: Deviation of surface temperature from its 
mean value as a function of latitude for same model as Fig. [3] 
Right panel: surface meridional flow. Bq = 200 G. 



The magnetic field structure and amplitude are illustrated 
by Figs. |5] and |6] The poloidal field amplitude, B = 200 G, 
is the mean value. Its maximum strength in the polar region is 
about 2.5 times larger. It is not poloidal field, however, which 
makes the principal influence on convective heat flux (jSJl and 
produces the temperature depletion of Fig. |4] The toroidal 
field of Fig.[6]is stronger by about one order of magnitude. 



Fig. 6. Toroidal magnetic field for same model as Fig. [3] 
Right panel shows radial profiles of the field for latitudes of 
20° (full line), 40° (dotted), and 60° (dashed). 



The anti-solar rotation of Fig.[3]can be interpreted along 
the following sequence. Toroidal magnetic field of Fig.[6]pro- 
duces thermal depletion of Fig.|4]by suppression of the con- 
vective heat flux 10. The thermal depletion drives a rela- 
tively fast meridional flow of Fig. via the barocline forc- 
ing of Eq. @. The meridional flow tends to make the angu- 
lar momentum uniform along its stream-lines in accord with 
the solution Q of the angular velocity equation (0) for the 
advection-dominated case thus producing relatively slow ro- 
tation of the equator. It should be noted, however, that con- 
vective fluxes of angular momentum (the A-effect, Rudiger 
& Hollerbach 2004 1 does also take part in formation of the 
differential rotation, and we return to the beginning of our 
interpretation by noticing that the toroidal magnetic field is 
produced from poloidal one by the differential rotation. 

The simulated anti-solar rotation of preceding section is 
the first theoretical indication for this type of differential rota- 
tion. Nevertheless, it is tempting to notice points of yet quali- 
tative agreement with observations of the giant star HD 3 1993 
by Strassmeier et al. ( 2003 1. The star shows differential rota- 
tion of about 12% of anti-solar type. The dark spots were 
observed on low latitudes and they had temperature contrast 
of about 200° K only which is compatible with (latitude- 
averaged) moderate temperature depletion of Fig.0] In a fur- 
ther agreement, the observed star had a warm pole. It is not 
clear, however, whether meridional flow was as fast as re- 
quired by our theoretical model to produce anti-solar rotation. 

Unfortunately, HD 31993 remains the only giant star for 
which the anti-solar type of rotation was detected incontro- 
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versially. The other candidate, HD 199178, is under debates 
(Hackmann et al. 120071 Petit et al. 2004). 

The anti-solar rotation was detected for six close binaries 
also (Strassmeier 2003). Consideration of section |5] suggests 
this type of star as another possibility for equatorial decelera- 
tion though no quantitative model for binaries was developed 
so far. 

5. What to observe? 

It seems to be a quite general theoretical statement that rel- 
atively fast meridional flow should be present whenever dif- 
ferential rotation is of anti-solar type. Accordingly, it is fun- 
damentally important to know whether observations can con- 
firm that all stars showing relatively slow equatorial rotation 
do simultaneously possess a large meridional flow. 

It is appropriate to estimate how fast the meridional cir- 
culation should be for producing the anti-solar differential 
rotation. Significance of meridional flow is defined by the 
Reynolds number, Re = u m R/v T ~ u^Rt/I 2 (£ and r are 
mixing length and time respectively). Judging from Fig. |2] 
the poleward meridional flow of 

u m > 30-^- (15) 
tR 

should be sufficient to support the equatorial deceleration. 
The fast circulation is expected to drive differential rota- 
tion of anti-solar type independently of whether the direction 
of the surface flow is towards pole or equator. Estimations 
by Riidiger ( 1989 1 suggest, however, that equatorward flow 
should be several times faster compared to Eq. dl5l to pro- 
duce anti-solar rotation. 

It may be noticed that the meridional circulation has been 
recognized as an important ingredient of stellar dynamos 
(Choudhuri et al. fl995l Dikpati & Gilman l200T1 Bonnano et 
al. l2002> . Accordingly, observational detections of the merid- 
ional flow would be helpful for stellar dynamo theory also. 

Another possibility for observational verification of the 
above theoretical predictions is to check whether all anti-solar 
rotators do indeed belong to one of two groups: spotted giants 
or close binaries. The suggestion that a fast meridional flow 
can be found only for the two groups of stars is, however, 
more speculative and less certain than the necessity of the 
meridional flow itself for rotation laws of anti-solar type. 

Measuring differential rotation and meridional flow on gi- 
ant stars can be a difficult observational task because of rela- 
tively long rotation periods of giants. The task may be a natu- 
ral subject for robotic astronomy because long observational 
series may be required. 

Acknowledgements. L.L.K. is grateful to A.I.P. for its hospitality 
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